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Abstract 



According to a formula that was put forward many decades ago 

the ground state energy per particle of an interacting, dilute Bose gas 

;h ■ at density p is 27rh'^pa/m to leading order in pa^ <C 1, where a is 

^ the scattering length of the interaction potential and m the particle 

mass. This result, which is important for the theoretical description of 
current experiments on Bose-Einstein condensation, has recently been 
established rigorously for the first time. We give here an account of 
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the proof that apphes to nonnegative, sphericahy symmetric potentials 
decreasing faster than 1/r^ at infinity. 

1 Introduction 

Recent progress in the trapping and cooling of atoms has made the ground 
state properties of dilute, interacting Bose gases accessible to experimental 
study [111 , [@] . In the theoretical description of such experiments an old for- 
mula for the ground state energy plays an important role. This formula, 
stated precisely in (|1.6|) below, is the subject of the present contribution, 
which is essentially an exposition of the paper , incorporating some new 
results from [Q and 0. 

We consider the Hamiltonian for A^ Bosons of mass m enclosed in a cu- 
bic box A of side length L and interacting by a spherically symmetric pair 
potential v{\xi — Xj\): 

N 
^7V = -/i^Ai+ Y, vi\xi-Xj\). (1.1) 

i=l i<i<j<N 

Here Xi & M.^, i = 1, . . . ,N are the positions of the particles, Aj the Lapla- 
cian with respect to Xi, and we have denoted fv^ /2m by /i for short. (By 
choosing suitable units /x could, of course, be eliminated, but we want to 
keep track of the dependence of the energy on Planck's constant and the 
mass.) The Hamiltonian (|1.1|) operates on symmetric wave functions in 
L'^{A^ ,dxi- ■ -dx^) as is appropriate for Bosons. The interaction poten- 
tial will be assumed to be nonnegative and to decrease faster than 1/r^ at 
infinity. 



We are interested in the ground state energy Eq{N,L) of (1J-) in the 
thermodynamic limit when A^ and L tend to infinity with the density p = 
N/L^ fixed. The energy per particle in this limit 

eo{p)= \imEo{pL^L)/{pL'). (1.2) 

Our results about eo(p) are based on estimates on Eo{N,L) for finite A^ 
and L, which are important, e.g., for the considerations of inhomogeneous 
systems in [Q. To define Eq{N,L) precisely one must specify the boundary 
conditions. These should not matter for the thermodynamic limit. To be 



on the safe side we use Neumann boundary conditions for the lower bound, 
and Dirichlet boundary conditions for the upper bound since these lead, 
respectively, to the lowest and the highest energies. 

For experiments with dilute gases the low density asymptotics of eo(p) is 
of importance. Low density means here that the mean interparticle distance, 
p-i/3 jg jjiuch larger than the scattering length a of the potential, defined as 

a = lim r - 1^, (1.3) 

where Uq solves the zero energy scattering equation, 

-2/iMo(r) + v{r)uo{r) = (1.4) 

with Uq{0) = 0. (The factor 2 in ( |1.4D comes from the reduced mass of the 
two particle problem.) Our main result is a rigorous proof of the formula 

eo(p) ~ 47r/ipa (1.5) 

for pa^ <^ 1, more precisely of 

Theorem 1.1 (Low density limit of the ground state energy) 

hm ^^ = 1. (1.6) 

pa3^0 ATTfXpa 

This formula is independent of the boundary conditions used for the definition 
of eo(p). 

The genesis of an understanding of eo(p) was the pioneering work [^] of 
Bogoliubov, and in the 50's and early 60's several derivations of ( |1.6| ) were 
presented [0], [§], even including higher order terms: 

^ = 1 + ^ipa^/' + 8 (f - v^) ipa') logipa') + 0{pa') (1.7) 

These early developments are reviewed in |^. They all rely on some special 
assumptions about the ground state that have never been proved, or on the 
selection of special terms from a perturbation series which likely diverges. 
The only rigorous estimates of this period were established by Dyson, who 
derived the following bounds in 1957 for a gas of hard spheres [|10| : 

1 ^ eo(p) ^ 1 + 2rV3 



10^2 " 47r/ipa " (1 - Fi/3)2 



with Y = Anpa^/S. While the upper bound has the asymptotically correct 
form, the lower bound is off the mark by a factor of about 1/14. But for 
about 40 years this was the best lower bound available! 

Since ( p. .61) is a basic result about the Bose gas it is clearly important 
to derive it rigorously and in reasonable generality, in particular for more 
general cases than hard spheres. The question immediately arises for which 
interaction potentials one may expect it to be true. A notable fact is that 
it not true for all v with a > 0, since there are two body potentials with 



positive scattering length that allow many body bound states ||Tl[. Our 
proof, presented in the sequel, works for nonnegative v, but we conjecture 
that (|1.6| ) holds if a > and v has no A^-body bound states for any A^. 
The lower bound is, of course, the hardest part, but the upper bound is not 
altogether trivial either. 

Before we start with the estimates a simple computation and some heuris- 
tics may be helpful to make ( |1.6| ) plausible and motivate the formal proofs. 

With uq the scattering solution and fo{r) = UQ{r)/r, partial integration 
gives 



/ {2fi\Vfo\' + v\mdx = An [ {2/xK(r) - K(r)/r)]2 + t;(r)|no(r)]2}rfr 

J\x\<R Jo 

= 87rfia\uo{R)\yR^ ^SiTfia for i? -> oo, (1.9) 



if Mq is normalized so that /o(-R) — ^ 1 as i? — >■ oo. Moreover, for positive inter- 
action potentials the scattering solution minimizes the quadratic form in ( |1.9|) 
for each R with uo{0) = and uo{R) fixed as boundary conditions. Hence the 
energy £^0(2, L) of two particles in a large box, i.e., L ^ a, is approximately 
SvTyua/L^. If the gas is sufficiently dilute it is not unreasonable to expect that 
the energy is essentially a sum of all such two particle contributions. Since 
there are A^(A^-l)/2 pairs, we are thus lead to Eo{N, L) ^ An fxaN {N -I) / L^ , 
which gives ( p..6| ) in the thermodynamic limit. 

This simple heuristics is far from a rigorous proof, however, especially for 
the lower bound. In fact, it is rather remarkable that the same asymptotic 
formula holds both for 'soft' interaction potentials, where perturbation the- 
ory can be expected to be a good approximation, and potentials like hard 
spheres where this is not so. In the former case the ground state is approx- 
imately the constant function and the energy is mostly potential: According 
to perturbation theory Eo{N, L) ^ A^(A^ — 1)/(2L'^) J v{\x\)dx. In particular 
it is independent of fi, i.e. of Planck's constant and mass. Since, however. 



J v{\x\)dx is the first Born approximation to Svr/xa (note that a depends on 
/i!), this is not in confiict with (|1.6| ). For 'hard' potentials on the other hand, 
the ground state is highly correlated, i.e., it is far from being a product of 
single particle states. The energy is here mostly kinetic, because the wave 
function is very small where the potential is large. These two quite different 
regimes, the potential energy dominated one and the kinetic energy dom- 
inated one, cannot be distinguished by the low density asymptotics of the 
energy. Whether they behave differently with respect to other phenomena, 
e.g., Bose-Einstein condensation, is not known at present. 

Bogolubov's analysis presupposes the existence of Bose-Einstein con- 
densation. Nevertheless, it is correct (for the energy) for the one-dimensional 
delta- function Bose gas [|I2|, despite the fact that there is (presumably) no 



condensation in that case. It turns out that BE condensation is not really 
needed in order to understand the energy. As we shall see, 'global' con- 
densation can be replaced by a 'local' condensation on boxes whose size is 
independent of L. It is this crucial understanding that enables us to prove 
Theorem 1.1 without having to decide about BE condensation. 

An important idea of Dyson was to transform the hard sphere potential 
into a soft potential at the cost of sacrificing the kinetic energy, i.e., effectively 
to move from one regime to the other. We shall make use of this idea in our 
proof of the lower bound below. But first we discuss the simpler upper bound, 
which relies on other ideas from Dyson's beautiful paper [^ . 



2 Upper bound 

The following generalization of Dyson's upper bound holds [Q, |^: 

Theorem 2.1 (Upper bound) Define pi = {N-l)/L^ andh = (47rpi/3)-^/^ 
For nonnegative potentials v, and b > a the ground state energy of ( \1 . 1\ ) with 
periodic boundary conditions satisfies 

Eo{N, L)/N < Aixppia ' /"^ „ ^ ^"^ . (2.1) 

(1-f) 

For Dirichlet boundary conditions the estimate holds with (const.) /L^ added 
to the right side. Thus in the thermodynamic limit and for all boundary 



conditions 



eoip) ^ l-Yy^ + Y^I^-\Y _ ^2^2^ 



47r/ipa - (1 - rV3) 

provided Y = 47rpa^/3 < 1. 



Remark. The bound ( p.l| ) holds for potentials with infinite range, provided 
b > a. For potentials of finite range Rq it can be improved for b > Rq to 

E,iN, L)/N < ATTf^p.a ^-^ ^^. (2.3) 

(1-f) 

Proof. We first remark that the expectation value of ( p..l| ) with any trial 
wave function gives an upper bound to the bosonic ground state energy, even 
if the trial function is not symmetric under permutations of the variables. 
The reason is that an absolute ground state of the elliptic differential operator 
( |1.1| ) (i.e. a ground state without symmetry requirement) is a nonnegative 
function which can be be symmetrized without changing the energy because 
(|1.1| ) is symmetric under permutations. In other words, the absolute ground 
state energy is the same as the bosonic ground state energy. 



Following ||I0[ we choose a trial function of the following form 



*(xi,...,XAf) = Fi(xi) ■ F2{xi,X2) ■■■FN{xi,...,Xn). (2.4) 

More specifically, Fi = 1 and Fj depends only on the distance of x, to 
its nearest neighbor among the the points Xi, . . . ,Xj_i (taking the periodic 
boundary into account): 

Fi{xi,...,Xi) = f{ti), ti=mm{\xi-Xj\,j = l,...,i-l) , (2.5) 

with a function / satisfying 

< / < 1, /' > 0. (2.6) 



The intuition behind the ansatz ( |2.4D is that the particles are inserted into the 
system one at the time, taking into account the particles previously inserted. 
While such a wave function cannot reproduce all correlations present in the 
true ground state, it turns out to capture the leading term in the energy 
for dilute gases. The form (|2.5|) is computationally easier to handle than an 

6 



ansatz of the type n^j/d^* ~ ^il)' which might appear more natural in 
view of the heuristic remarks at the end of the last section. 
The function / is chosen to be 

f^^^^,foir)/Ub) forO<r<6, ^^.7) 

1 for r > b, 



with fo{r) = uo{r)/r. The estimates (p.l|) and ( p.3|) are obtained by some 



what lengthy computations similar as in |Ty], but making use of ( |1.9| ). For 
details we refer to [Q and [^. 

A test wave function with Dirichlet boundary condition may be obtained 
by localizing the wave function (|2.4|) on the length scale L. The energy cost 
per particle for this is (const. )/L^. D 

3 Lower bound 

To get an idea why the lower bound for the bosonic ground state energy of 
(|1 . 1| ) is not easy to obtain let us consider the relevant length scales of the 
problem. These are 

• The scattering length a. 

• The mean particle distance p~^/^. 

• The 'uncertainty principle length' ic, defined by /i£~^ = eo(p), i.e., 
4 ~ (pa)-'/'. 

The length ic is sometimes called 'correlation length' or 'healing length'. 
The name 'uncertainty principle length' is justified by the fact that this is 
the shortest length scale on which the bosons can be localized without raising 
the energy per particle above cq, according to the uncertainty principle. For 
dilute gases pa^ -C 1 and hence 

a « p~i/3 ^ ^pa^yy^p-y^ ^ 4. (3.1) 

Bosons in their ground state are therefore 'smeared out' over distances large 
compared to the mean particle distance and their individuality is entirely 
lost. Fermions, on the other hand, prefer to sit in private rooms, i.e., ic can 
be comparable to p~^^^. In this respect the quantum nature of Bosons is 



much more pronounced than for Fermions. The three different length scales 
for Bosons will play a role in the proof below. 
Our lower bound for eo(p) is as follows. 

Theorem 3.1 (Lower bound in the thermodynamic limit) For a pos- 
itive potential v with finite range and Y small enough 

"""^^^ >{1-CY'/'') (3.2) 



Anupa 



with C a constant. If v does not have finite range, hut decreases at least as 
fast as Ijr^^^ at infinity with some e > 0, then an analogous hound to ^3.2J 
holds, hut with C replaced hy another constant and 1/1 7 hy another exponent, 
hoth of which may depend on e. 



It should be noted right away that the error term —CY^^^'^ in (|3.2| ) is of 



no fundamental significance and is not believed to reflect the true state of 
affairs. Presumably, it does not even have the right sign. We mention in 
passing that C can be taken to be 8.9 0. 

As mentioned in the Introduction a lower bound on Eq{N, L) for finite N 
and L is of importance for applications to inhomogeneous gases, and in fact 
we derive (|3.2| ) from such a bound. We state it in the following way: 



Theorem 3.2 (Lower bound in a finite box) For a positive potential v 
with finite range there is a 6 > such that the the ground state energy of 
( \1.1\) with Neumann conditions satisfies 

Eo{N,L)/N>A'Kfipa{l-CY^/^'^) (3.3) 

for all N and L with Y < 5 and L/a > C'Y^^I^'^ . Here C and C are 
constants, independent of N and L. (Note that the condition on L/a requires 
in particular that N must he large enough, N > {coiast.)Y~^'^'^ .) As in 
Theorem |5[J such a hound, hut possihly with other constants and another 
exponent for Y , holds also for potentials v of infinite range decreasing faster 
than 1/r^ at infinity. 



The first step in the proof of (|3.2| ) is a generalization of a lemma of Dyson, 
which allows us to replace f by a 'soft' potential, at the cost of sacrificing 
kinetic energy and increasing the effective range. 



Lemma 3.3 Letvir) > with finite range Rq. Let U{r) > be any function 
satisfying J U{r)r'^dr < 1 and U{r) = for r < Rq. Let B C R^ be star 
shaped with respect to (e.g. convex with G B). Then for all differentiable 
functions ip 

'"[/i|VV^P + it;|V^P]>/ia/"[/|^p. (3.4) 

B Jb 



Proof. Actually, ( p. 41) holds with yu|V0(x)p replaced by the (smaller) radial 
kinetic energy, fj,\d(f){x) / dr]"^ , and it suffices to prove the analog of ( p.4|) for 
the integral along each radial line with fixed angular variables. Along such a 
line we write 0(x) = u{r)/r with m(0) = 0. We consider ffist the special case 
when when [/ is a delta-function at some radius R > Rq, i.e., 

U(r) = ^5ir ~ R). (3.5) 



For such U the analog of (3]^) along the radial line is 



\^i[u'{r)~{u{r)/r)]' + \v{r)\u{r)]'}dr>r , ^^., . ^, '[oV^zf 
\^a\u{R)\^/R'' It R < Ri 

(3.6) 

where Ri is the length of the radial line segment in B. The case Ri < R is 
trivial, because fi\dip/dr\'^ + |v|V'P > 0. (Note that positivity of v is used 
here.) li R < Ri we consider the integral on the the left side of ( p.6|) from 
to -R instead of -Ri and minimize it under the boundary condition that 
m(0) = and u{R) is a fixed constant. Since everything is homogeneous in u 
we may normalize this value to u{R) = R — a. This minimization problem 



leads to the zero energy scattering equation (p..4|). Since v is positive, the 
solution is a true minimum and not just a stationary point. 

Because f (r) = for r > _Ro the solution, uq, satisfies Uo{r) = r — a for 
r > Rq. By partial integration, 

R 

{^K(r) - (Mo(r)/r)]2 + |t;(r)|Mo(r)]2}rfr = /ia|i? - a|V^'- (3.7) 



But l-R — a\^ / R^ is precisely the right side of ( p.6|) if u satisfies the normal- 
ization condition. 



This derivation of ( ^.4| ) for the special case ( |3.5| ) imphes the general case, 
because every U can be written as a superposition of (5-functions, U{r) = 
J R-^6{r - R) U{R)R^dR, and / U{R)R^dR < 1 by assumption. D 

By dividing A for given points xi, . . . , xn into Voronoi cells Bi that contain 
all points closer to Xi than to Xj with j y^ i (these cells are star shaped w.r.t. 



Xi, indeed convex), the following corollary of Lemma 3^ can be derived in 
the same way as the corresponding Eq. (28) in [|10|] . 



Corollary 3.4 For any U as in Lemma [J.Jl 

Hn > liaW (3.8) 

with 



N 



iy(fi,...,x^) = ^?7(t,), (3.9) 

where ti is the distance of Xi to its nearest neighbor among the other points 
Xj, j = 1,...,N, I.e., 

ti{xi, . . . ,Xj\[) = min \xi — Xj\. (3.10) 



(Note that ti has here a slightly different meaning than in ( ^.51 ), where it 
denoted the distance to the nearest neighbor among the Xj with j < i — 1.) 
Dyson considers in []TU[ a one parameter family of f/'s that is essentially 



the same as the following choice, which is convenient for the present purpose: 

otherwise. 



We denote the corresponding interaction ( |3.9| ) by Wr. For the hard core gas 
one obtains 

E{N, L) > sup inf fiaWnixi, ...,xn) (3.12) 

where the infimum is over {xi, . . . , x^) G A^ with \xi — Xj\ > Rq = a, because 
of the hard core. At fixed R simple geometry gives 

inf^ Wr{x^,...,xn)> f4--|^l (3-13) 

ixi,...,XN) yR"^ pR J 

10 



with certain constants A and B. An evaluation of these constants gives 
Dyson's bound 

E{N,L)/N>^A7ifxpa. (3.14) 

The main reason this method does not give a better bound is that R 
must be chosen quite big, namely of the order of the mean particle distance 
p~^/^, in order to guarantee that the spheres of radius R around the N points 
overlap. Otherwise the infimum of Wr will be zero. But large R means that 
Wr is small. It should also be noted that this method does not work for 
potentials other than hard spheres: If \xi — Xj\ is allowed to be less than Rq, 
then the right side of (|3.12| ) is zero because U{r) = for r < Rq. 



For these reasons we take another route. We still use Lemma |3.4| to get 



into the soft potential regime, but we do not sacrifice all the kinetic energy 
as in (13.8). Instead we write, for e > 



-') 



Hj^ = eHn + (1 - €)Hn > sTn + (1 - e)HN (3.15) 

with T/v = — Ylii ^i ^^^ US6 (|378|) only for the part (1 — e)H^. This gives 

Hn > sTn + (1 - e)iiaWB. (3.16) 

We consider the operator on the right side from the viewpoint of first order 
perturbation theory, with eT/v as the unperturbed part, denoted Hq. 

The ground state of Hq in a box of side length L is '^o{xi, . . . , xn) = 
]^-3N/2 ^^^ ^g denote expectation values in this state by (•)o. A computation, 
cf. Eq. (21) in p, gives 

47rp(l-^) > {Wn)o/N (3.17) 

> 4vrp(l-i)(l-f)^(l + 4vrp(l-^)(/?3_/23)/3))-i, 

The rationale behind the various factors is as follows: (1 — -^) comes from the 
fact that the number of pairs is A^(iV- 1)/2 and not N'^/2, (1 - 2R/Lf takes 
into account the fact that the particles do not interact beyond the boundary 
of A, and the last factor measures the probability to find another particle 
within the interaction range of the potential Ur for a given particle. 

The first order result ( 3.17| ) looks at first sight quite promising, for if we 



let L — > oo, A^ — *> oo with p = N/L fixed, and subsequently take i? — *■ co, 

11 



then {Wr)qIN converges to 47rp, which is just what is desired. But the 



first order resuh ( p.l7|) is not a rigorous bound on Eq{N,L), we need error 
estimates, and these will depend on e, R and L. 

We now recall Temple 's inequality [|^ for the expectations values of an 
operator H = Hq+V in the ground state {■)o of Hq. It is a simple consequence 
of the operator inequality 

{H -Eo){H -Ei)>0 (3.18) 

for the two lowest eigenvalues, Eq < Ei, oi H and reads 

E«>{H),-!^^ghzAm (3.19) 

provided Ei — {H)q > 0. Furthermore, if ^ > we may use Ei > E^ — 
second lowest eigenvalue of Hq and replace Ei in (|3.19| ) by _EJ . 
From ( p. 171 ) and ( |3.19| ) we get the estimate 



Eo{N,L) 
with 



A^ 



>ATrfiap{l-£{p,L,R,e)) (3.20) 



l-£{p,L,R,e) = (l-e){l-^){l-fY{l + fp{l-j^)iR'-Rl))) 

{Wn)o{Ei''^ - fia{WR)o) 
To evaluate this further one may use the estimates (|3.17| ) and the bound 



(3.21) 



{W'n)o<S^^{Wn)o (3.22) 

which follows from U"^ = 3{R^ — Rq)^^Ur together with the Cauchy-Schwarz 
inequality. A glance at the form of the error term reveals, however, that it 
is not possible here to take the thermodynamic limit L ^ oo with p fixed: 
We have E[ = evr/i/L^ (this is the kinetic energy of a single particle in 
the first excited state in the box), and the factor El — p.a{WR)o in the 
denominator in ( p.21| ) is, up to unimportant constants and lower order terms. 



~ {eL ^ — ap'^L^). Hence the denominator eventually becomes negative and 
Temple's inequality looses its validity if L is large enough. 

12 



As a way out of this dilemma we divide the big box A into cubic cells of 
side length i that is kept fixed as L — > oo. The number of cells, L^/i^, on the 
other hand, increases with L. The N particles are distributed among these 
cells, and we use ( |3.21|) , with L replaced by i, N by the particle number, n, 
in a cell and p by n/i^, to estimate the energy in each cell with Neumann 
conditions on the boundary. This boundary condition leads to lower energy 
than any other boundary condition. For each distribution of the particles we 
add the contributions from the cells, neglecting interactions across bound- 
aries. Since w > by assumption, this can only lower the energy. Finally, 
we minimize over all possible choices of the particle numbers for the various 
cells adding up to A^. The energy obtained in this way is a lower bound to 
Eq{N, L), because we are effectively allowing discontinuous test functions for 
the quadratic form given by H^- 

In mathematical terms, the cell method leads to 

Eo{N,L)/N > {pe^)-HniJ2cnEo{n,i) (3.23) 

n>0 

where the infimum is over all choices of coefficients c„ > (relative number 
of cells containing exactly n particles), satisfying the constraints 

5^ c„ = 1, Yl C"^ = P^^- (3-24) 

n>0 n>0 

The minimization problem for the distributions of the particles among 
the cells would be easy if we knew that the ground state energy EQ{n, i) (or 
a good lower bound to it) were convex in n. Then we could immediately 
conclude that it is best to have the particles as evenly distributed among the 
boxes as possible, i.e., Cn would be zero except for the n equal to the integer 
closest to pi^. This would give 

E^i^ > 4rrpap (1 - S{p, i, R, e)) (3.25) 

i.e., replacement of L in ( p. 201) by i, which is independent of L. The blow up 
of S for L ^ oo would thus be avoided. 

Since convexity of Eo{n,i) is not known (except in the thermodynamic 
limit) we must resort to other means to show that n = 0{pi^) in all boxes. 
The rescue comes from superadditivity of Eo{n,i), i.e., the property 

Eo{n + n\ i) > Eoin, i) + Eo{n' , £) (3.26) 

13 



which follows immediately from f > by dropping the interactions between 
the n particles and the n' particles. The bound ( p.26| ) implies in particular 
that for any n,p eN with n > p 



n 



E{nJ)>[n/p]E{pJ)>—E{p, 



(3.27) 



since the largest integer [n/p] smaller than n/p is in any case > n/{2p). 

The way ( |3.27D is used is as follows: Replacing L hj i, N hj n and p by 
n/i^ in ( |3.2(]| ) we have for fixed R and e 



EoinJ)>^!^n{n-l)Kin,i) 



(3.28) 



with a certain function K(n, i) determined by ( ^.211 ). We shall see that K is 
monotonously decreasing in n, so that if p G N and n < p then 



Eo{n,£) > ^n{n - l)K{p,i). 



(3.29) 



We now split the sum in ( p.23| ) into two parts. For n < p we use ( p.29|) , and 
for n > p we use (|3.27|) together with ( p.29|) for n = p. The task is thus to 
minimize 



y^ Cnn{n - 1) + I ^ c„n(p - 1) 

n>p 



(3.30) 



n<p 



subject to the constraints (|3.24|) . Putting 

k := pi^ and t := >^ CnU < k 



(3.31) 



n<p 



we have '^^^^p CnU = k — t, and since n(n — l) is convex in n, and J2n<p ^n < 1 



the expression (|3.30|) is 



>t(t-i) + Uk-t){p-i). 



(3.32) 



We have to minimize this for 1 < t < A;, li p > 4k the minimum is taken at 
t = k and is equal to k{k — 1). Altogether we have thus shown that 



Eo{N,L] 
N 



> Anpap M - ^ j K{Apf 



(3.33) 
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What remains is to take a closer look at K{4pi'^,i), which depends on 
the parameters e and R besides i, and choose the parameters in an optimal 
way. /^From ( |3.21| ) and ( p.22| ) we obtain 



K{n,i) = (l-e)(l-f)'(l + fp(l-i)(i?^-i?o')))"' 

1_3 an \ 

TT (i?3 _ Rl){ei-^ - Aai-^n -1))/ 

The estimate (|3.28|) with this K is valid as long as the denominator in the 
last factor in ( ^.34} ) is > 0, and in order to have a formula for all n we can 
take as a trivial lower bound in other cases or when (|3.28|) is negative. 



As required for ( |3.29| ), K is monotonously decreasing in n. We now insert 
n = 4pi^ and obtain 



K{Apf,i) > {l-E){l-^f {1 + {const. )Y{i/aYiR^ -Rl)/i^) ' 

i^ (const. )F 



(i?3 _ Rl) {e{a/iy - (const. )r2(£/a)3) J (^.35) 

with Y = 4:7ipa^/3 as before. Also, the factor 

(l - ^) = (1 - {const.)Y~\a/if) (3.36) 

in ( ^.33| ) (which is the ratio between n(n — 1) and n'^) must not be be for- 
gotten. We now make the ansatz 

e~F", a/i^Y^, {R^-Rl)/f ^Y^ (3.37) 

with exponents a, (3 and 7 that we choose in an optimal way. The conditions 
to be met are as follows: 

• e{a/(iY — (const.)F^(£/a)3 > 0. This holds for all small enough F, 
provided a + 5/? < 2 which follows from the conditions below. 

• a > in order that e — > for y — > 0. 

• 3/3 - 1 > in order that Y-^{a/tj^ -^ for for Y ^ 0. 

• l-3/3 + 7>0in order that Y{i/aY{R^ - Rl)/i^ -> for for Y ^ 0. 

• 1 — a — 2/9 — 7>0to control the last factor in ( 3.35|) . 
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Taking 

a = 1/17, P = 6/17, 7 = 3/17 (3.38) 

all these conditions are satisfied, and 

a = 3/3-l = l-3/3 + 7 = l-a-2/5-7 = 1/17. (3.39) 

It is also clear that 2R/i ~ Y^^^ = F^/^^, up to higher order terms. This 
completes the proof of Theorems 3.1 and 3.2, for the case of potentials with 
finite range. By optimizing the proportionality constants in ( 3.37|) one can 



show that C = 8.9 is possible in Theorem 1.1 |^. The extension to potentials 
of infinite range decreasing faster than 1/r^ at infinity is obtained by approx- 
imation by finite range potentials, controlling the change of the scattering 
length as the cut-off is removed. See Appendix B in |^ for details. A slower 
decrease than 1/r^ implies infinite scattering length. D 

The exponents ( p.38| ) mean in particular that 

a<^R<^ p-^/^ < £ < (pa)-^/^, (3.40) 

whereas Dyson's method required R ~ p~^^^ as already explained. The 
condition p~^/^ <C ^ is required in order to have many particles in each box 
and thus n{n — 1) Ri n'^. The condition i <^ (pa)~^^^ is necessary for a 
spectral gap gap ^ eo(p) in Temple's inequality. It is also clear that this 
choice of £ would lead to a far too big energy and no bound for eo(p) if we 
had chosen Dirichlet instead of Neumann boundary conditions for the cells. 
But with the latter the method works! 
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